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^ ! Abstract 

We calculate the binding energies for multi-nucleon bound states with the nuclear mass number 
less than or equal to 4 in 2+1 flavor QCD at the lattice spacing of a = 0.09 fm employing a 
relatively heavy quark mass corresponding to = 0.51 GeV. To distinguish a bound state from 
attractive scattering states, we investigate the volume dependence of the energy shift between the 
ground state and the state of free nucleons by changing the spatial extent of the lattice from 2.9 fm 
to 5.8 fm. We conclude that ^He, ^He, deuteron and dineutron are bound at = 0.51 GeV. We 
compare their binding energies with those in our quenched studies and also with several previous 
investigations. 

PACS numbers: 11.15.Ha, 12.38.Aw, 12.38.-t 12.38.Gc 
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I. INTRODUCTION 



Lattice QCD has a potential ability to quantitative understand the nature of nuclei, 
whose characteristic feature is a hierarchical structure in the strong interaction. The nuclear 
binding energy is experimentally known to be about 10 MeV per nucleon, which is much 
smaller than the typical energy scale of hadrons. A measurement of the binding energies is 
therefore the first step for direct investigation of nuclei in lattice QCD. A key ingredient in 
the study is a systematic change of the spatial volume of the lattice to distinguish a bound 
state from an attractive scattering state. 

We carried out a first attempt to measure the binding energies of the ^He and ^He nuclei 
in quenched QCD with a rather heavy quark mass corresponding to m^r = 0.80 GeV, thereby 
avoiding heavy computational cost l|. We followed this work by a renewed investigation 
of the bound state for the two-nucleon channel in quenched QCD at the same quark mass, 
which found that not only the deuteron in the ^Si channel but also the dineutron in the ^Sq 
channel is bound Independently, NPLQCD collaboration reported the possibility that a 
bound state is formed in both channels at m^r = 0.39 GeV in 2+1 flavor QCD [3|]. They later 
confirmed the bound states for the helium nuclei and the two-nucleon channels at = 0.81 
GeV in 3-fiavor QCD taking a different choice for the quark and gluon actions j4|. 

In this article we report on our investigation of the dynamical quark effects on the binding 
energies of the helium nuclei, the deuteron and the dineutron. We perform 2+1 flavor lattice 
QCD simulation with the degenerate up and down quark mass corresponding to m^r = 0.51 
GeV. Four lattice sizes are employed to take the infinite spatial volume limit: 32^ x 48, 
40^ X 48, 48^ X 48and 64^ x 64, whose spatial extent ranges from 2.9 fm to 5.8 fm with the 
lattice spacing of a = 0.08995(40) fm Q. 

For the helium nuclei our main interest lies in the magnitude of the binding energies, 
since all studies carried out so far, both in quenched and in unquenched QCD and for 
several quark mass values, agree on the bound state nature for helium nuclei. Much more 
intriguing is the two-nucleon system, for which there are two ways to study. One is a direct 
investigation j2-4, 6-9| in which one calculates the two-nucleon Green's functions directly in 



lattice QCD, and the other is an indirect calculation by means of the two-nua 



potential extracted from the two-nucleon wave function in lattice QCD 10|, lll |. 



eon effective 



So far only the former method has reported the binding energies of the two-nucleon 
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systems. In quenched QCD the bound state nature has been confirmed for both channels 
at = 0.80 GeV in our recent work On the other hand, unquenched studies show a 
comphcated situation. A somewhat early study in 2+1 flavor QCD with a mixed action js] 
reported a positive energy shift (repulsive interaction) in both channels at ttItt < 0.59 GeV. 
More recently, however, deep bound states were observed at = 0.81 GeV in 3-flavor 
QCD We hope to shed light on this situation with our own investigation in 2+1 flavor 
QCD. 

This paper is organized as follows. In Sec. [Ill we explain simulation details including 
the simulation parameters and the interpolating operators for the multi-nucleon channels. 
Section IIIII presents the results of the binding energies for the helium nuclei, the deuteron 
and the dineutron. We compare our results with those in the previous studies. Conclusions 
and discussions are summarized in Sec. IIVI 



II. SIMULATION DETAILS 



A. Simulation parameters 



We generate 2+1 flavor gauge configurations with the Iwasaki gauge action 12[ and the 
non-perturbative 0(a)-improved Wilson quark action at /3 = 1.90 with csw = 1-715 13|. 
The lattice spacing is a = 0.8995(40) fm, corresponding to a^^ = 2.194(10) GeV, determined 
with mn = 1.6725 GeV Q. We take four lattice sizes, x T = 32^ x 48, 40^ x 48, 48^ x 48 
and 64^ x 64, to investigate the spatial volume dependence of the ground state energy shift 
between the multi-nucleon system and the free nucleons. The physical spatial extents are 
2.9, 3.6, 4.3 and 5.8 fm, respectively. Since it becomes harder to obtain a good signal-to- 

nn 

noise ratio at lighter quark masses for multi-nucleon systems p, HJ], we employ the hopping 
parameters {nud, i^s) = (0.1373316,0.1367526) which correspond to m.,r = 0.51 GeV and 
m^v = 1.32 GeV and the physical value for the strange quark mass. These values are chosen 



based on the previous results for and obtained by the PACS-CS Collaboration J5 
We employ the domain-decomposed Hybrid-Monte- Carlo (DDHMC) algorithm Il6 
for the degenerate light quarks and the UV- filtered PHMC (UVPHMC) algorithm |l8 



the strange quark employing the Omelyan-Mryglod-Folk integrator 19 



20[ . The algorithmic 



details are given in Ref. 15l|. We summarize simulation parameters in Table H] including the 
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block sizes in DDHMC and the polynomial order in UVPHMC. We take r = 1 for the 
trajectory length of the molecular dynamics in all the runs. The step sizes are chosen such 
that we obtain reasonable acceptance rates presented in Table [H We generate the gauge 
configurations in a single run except for the L = 64 case for which we carry out two runs. 
The total trajectory length is about 2000 for all the volumes, except 4000 for the case of 
smallest volume. 

B. Calculation method 

We extract the ground state energies of the multi-nucleon systems and the nucleon state 
from the correlation functions 



with O being appropriate operators for ^He, '^He, two-nucleon ^Si and ^Sq channels, and the 
nucleon state (see the next subsection for actual expressions). 

We carry out successive measurements in the interval of 10 trajectories. The errors are 
estimated by jackknife analysis choosing 200 trajectories for the bin size for all volumes, 
except for the largest volume for which we use 190. The numbers of configurations are listed 
in Table Hn We attempt to extract as much information as possible from each configuration by 
repeating the measurement of the correlation functions for a number of sources at different 
spatial points and time slices. For the 48^ and 64^ lattices, we calculate the correlation 
functions not only in the temporal direction but also in the three spatial directions exploiting 
the space-time rotational symmetry. We found that this procedure effectively increases 
statistics by a factor of four. This factor is included in the number of measurements on each 
configuration given in Table [TTl 

We are interested in the energy shift between the ground state of the multi-nucleon system 
and the free nucleons on an box. 



with Eq being the lowest energy level for the multi-nucleon channel, the number of 
nucleon and m^v the nucleon mass. This quantity is directly extracted from the ratio of the 
multi-nucleon correlation function divided by the A^Ar-th power of the nucleon correlation 



Go{t) = {0\O{t)Om0) 



(1) 



AEl = Eo- NNrriN 



(2) 



4 



function 

where the same source operator is chosen for the numerator and the denominator. We also 
define the effective energy shift as 

AEf = \n( -^^] , (4) 

which is useful to check the plateau region in later sections. 

Note that the definitions for AEl and AEf^ follow those in Ref. but are opposite to 
those used in Ref. 



C. Interpolating operators 

We use an interpolating operator for the proton given by 

Pa = £abci[UaYC-f5db)u^, (5) 

where C = 7472 and a and a, b, c are the Dirac index and the color indices, respectively. 
The neutron operator Ua is obtained by replacing u° by (i° in the proton operator. To save 
the computational cost we use the non-relativistic quark operator, in which the Dirac index 
is restricted to the upper two components. 

The ^He nucleus has zero total angular momentum, positive parity = O"*" and zero 
isospin 1 = 0. We employ the simplest "^He interpolating operator with zero orbital angular 
momentum L = 0, and hence J = S with S being the total spin. Such an operator was 



already given long time ago in Ref. 21 1 



^He = (xr/ - xv) , (6) 



where 



X = ^([+ - +-] + [- + -+] - [+ - -+] - [- + +-]), (7) 
X = - +-] + [- + -+] + [+ - -+] + [- + +-] - 2[+ + — ] - 2[- - ++])(8) 



with +/— being up/down spin of each nucleon, and 77, r/ are obtained by replacing +/— in 
X,X by p/n for the isospin. Each nucleon in the sink operator is projected to zero spatial 
momentum. 
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We also calculate the correlation function of the ^He nucleus whose quantum numbers 



are = , I = ^ and = |. We employ the interpolating operator in Ref. 22 1 



^He = {\p.n+p+) - \p+n+p^) + \n+p+p^) - \n+p_p+) + \p+p^n+) - \p_p+n+)) , (9) 

with the zero momentum projection on each nucleon in the sink operator. 
The two-nucleon operators for the ^Si and ^Sq channels are given by 

NN^sAt) = ^ [p+{tM) - n+(tK(t)] , (10) 

iViViSo(t) = ^ [n4t)n-{t) - n-it)n+{t)] . (11) 

In the spin triplet channel the operators for the other two spin components are constructed 
in a similar way. We take average over the three spin components. 

The quark propagators are solved with the periodic boundary condition in all the spatial 
and temporal directions using the exponentially smeared source of form 

q'{x,t) = J2^e~^^'~Hy,t) (12) 
y 

after the Coulomb gauge fixing. We choose the smearing parameters depending on the 
volume (see Table [Tll) in order to obtain reasonable plateaux of the effective energy for 
the ground states in the multi-nucleon channels as well as for the nucleon. For the source 
operators explained above we insert the smeared quark fields of Eq. f lT^ for each nucleon 
operator located at the same spatial point x. Each nucleon in the sink operator, on the 
other hand, is composed of the point quark fields, and projected to zero spatial momentum. 

D. Difficulties for multi-nucleon channel 

There are several computational difficulties in the calculation of the correlation functions 
Goit) for the ^He and ^He channels. One is a factorially large number of Wick contractions 
for the quark- ant iquark fields. A naive counting gives {2Np + Nn)\{2Nn + Np)\ for a nucleus 
composed of Np protons and Nn neutrons, which quickly becomes prohibitively large beyond 
three-nucleon systems, e.g., 2,880 for '^He and 518,400 for ^He. To overcome the difficulty. 



we use the reduction techniques proposed in our exploratory work [l|. After the reduction, 
only 1107 (93) contractions are required for the correlation function in the '^He (^He) channel. 
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Other reduction techniques for the large number of the Wick contractions have been proposed 
for the multi-meson 23| and multi-baryon 2J, |25| channels. 



Another difficulty in studying a multi-nucleon bound state is the identification of the 
bound state nature in a finite volume, because an attractive scattering state yields a similar 
energy shift due to the finite volume effect 



26 



28 



. To solve the problem we need to inves- 
tigate the volume dependence of the measured energy shift Q]: For a scattering state, 
the energy shift decreases in proportion to l/L"^ at the leading order in the 1/L expan- 



sion 



26 



29| , while for a bound state the energy shift remains at a finite value in the infinite 
spatial volume limit. In order to distinguish a non-zero constant from a behavior in 

the energy shift, we employ four spatial extents from 2.9 to 5.8 fm. 



III. RESULTS 
A. Nucleon 

We first show the effective nucleon mass on the (5.8 fm)'^ box in Fig. [1] as a typical result. 
The plateau of the effective mass is clearly observed. A fit result of the correlation function 
with an exponential form is also drawn in the figure with the one standard deviation error 
band. We list the nucleon mass together with the pion mass in Table HIl 



B. ^He nucleus 

The effective energy shift Ai?£^ defined in Eq. (jlj) is plotted in Fig. [21 The signal is 
clear up to t = 12, beyond which the statistical error increases rapidly. The energy shift 
AEi is extracted from R{t) of Eq. ([3]) by an exponential fit over the range oi t = 10-14. 
The fit result is denoted by the solid lines with the statistical error band in Fig. [2l The 
systematic error in the fit is estimated from the variation of the fit results with the minimum 
or maximum time slice changed by ±1. Results with similar quality are obtained on other 
volumes. We summarize the values of AE'l with the statistical and systematic errors in 
Table [ml 

Figure [3l shows the volume dependence of Ai?^ as a function of The inner bar 

denotes the statistical error and the outer bar represents the statistical and systematic errors 
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combined in quadrature. The negative energy shifts are obtained in all the four volumes. 
We extrapolate the results to the infinite volume limit with a simple linear function of l/L'^, 

AEl = AE^ + ^. (13) 

The systematic error is estimated from the variation of the results obtained by alternative 
fits which contains a constant fit of the data and a fit of the data obtained with a different fit 
range in t. The non-zero negative value obtained for the infinite volume limit AE^q shown 
in Fig. [3] and Table IIIII leads us to conclude that the ground state is bound in this channel 
for the quark masses employed. The binding energy —AEoo = 43(12) (8) MeV, where the 
first error is statistical and the second one is systematic, is consistent with the experimental 
result of 28.3 MeV and also with the previous quenched result at = 0.80 GeV Note 
that the error is still quite large. 

A recent work in 3-fiavor QCD at 777,7^ = 0.81 GeV reported a value 110(20) (15) MeV 
for the binding energy of ^He nucleus j4|. This is about three times deeper than our value. 
Whether this difference can be attributed to the quark mass dependence in unquenched 
calculations needs to be clarified in future. 



C. ^He nucleus 

Figure m shows the effective energy shift AEf^ of Eq. (j4]). The quality of the signal is 
better than the ^He channel in Fig. |2l An exponential fit of R{t) in Eq. ([3]) with the range 
of t = 9-14 yields a negative value, which is denoted by the solid lines with the statistical 
error band in Fig. SI The systematic error in the fit is estimated in the same way as in the 
^He case. 

As listed in Table IIIII we find non-zero negative values for the energy shift AEl for all 
the volumes. The volume dependence is illustrated in Fig. |5] as a function of with the 
inner and outer error bars as explained in the previous subsection. We carry out a linear 
extrapolation of Eq. f lT5]) . The systematic error is estimated in the same way as in the ^He 
channel. The energy shift extrapolated to the infinite spatial volume limit is non-zero and 
negative, see Fig. [5] and Table llllt which means that the ground state is a bound state in 
this channel. The value of —AE^ = 20.3(4.0) (2.0) MeV is roughly three times larger than 
the experimental result, 7.72 MeV, though consistent with our previous quenched result at 
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= 0.80 GeV 



In 3-flavor QCD — Ai^oo = 71 (6) (5) MeV was reported [J] at a heavier quark mass 
corresponding to = 0.81 GeV. Here again future work is needed to see if a quark mass 
dependence explains the difference from the experiment. 

D. Two-nucleon channels 

1. Present work 

In Fig. |6] we show the time dependence for AEf^ of Eq. (jlj) in the '^Si channel. The 
signals are lost beyond t ^ 14. We observe negative values beyond the error bars in the 
plateau region of t = 9-14. We extract the value of AEl from an exponential fit for R{t) of 
Eq. (I3l) in the range of t = 9-14. The systematic error of the fit is estimated as explained in 
the previous subsections. 

Figure [7] shows the result for AE'j^ in the ^Sq channel on the (5.8 fm)'^ box. The value 
of AEf^ is again negative beyond the error bars in the plateau region, though the absolute 
value is smaller than in the ^Si case. The energy shift AE^ is obtained in the same way as 
for the ^Si channel. 

The volume dependences of AEl in the ^Si and ^Sq channels are plotted as a function of 
in Figs. M and |9l respectively. The numerical values of AEl on all the spatial volumes 
are summarized in Table llVj where the statistical and systematic errors are given in the first 
and second parentheses, respectively. There is little volume dependence for AEl, indicating 
a non-zero negative value in the infinite volume and a bound state, rather than the 
dependence expected for a scattering state, for the ground state for both channels. 

The binding energies in the infinite spatial volume limit in Table HV] are obtained by 
fittin g th e data with a function including a finite volume effect on the two-particle bound 



state 
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A^.^-^^ + f?>; - (14) 

where 7 and Cy are free parameters, n is three-dimensional integer vector, and Yl'n denotes 
the summation without |n| = 0. The binding energy —AE^o is determined from 

- AEoo = — , (15) 
rriN 



where we assume 

2\/^^-7'-2m;v~- — . (16) 

The systematic error is estimated from the variation of the fit results choosing different fit 
ranges in the determination of AEl and also using constant and linear fits as an alternative 
fit forms. We obtain the binding energies -AEoo = ll-5(l.l)(0.6) MeV and 7.4(1. 3)(0.6) 
MeV for the ^Si and ^Sq channels, respectively. The result for the "^Si channel is roughly 
five times larger than the experimental value, 2.22 MeV. Our finding of a bound state in the 
^Sq channel contradicts the experimental observation. These features are consistent with 
our quenched results with a heavy quark mass corresponding to m.^ = 0.80 GeV 



2. Comparison with previous studies 



A number of studies have been performed for the two-nucleon channel after the first work 
of Ref. [3]. It is therefore instructive to summarize the results and make a comparison with 
each other. Table IVl tabulates, in chronological order, the results for —AE^ for the ^Si and 
^Sq channels together with the pion mass and the spatial extent L in physical units. 
The numbers are plotted in Figs. [10] and [11] for the '^Si and ^Sq channels, respectively, as a 
function of m^. 

The early studies in Refs. fly , [ill employed a single volume, and we do not observe 
a common feature or trend among them. The positive values for —AEl in Ref. ^ means 
repulsive interaction for both channels, which is not seen in other studies. The results for 



-AEl in Ref. 



ll| is an order of magnitude smaller compared to other groups, probably due 



to significant contamination from excited states. 



The four recent studies 



have made a systematic investigation of the spatial volume 
dependence. Our quenched and 2+1 flavor results show qualitatively the same feature that 
the binding energy for the ^Si channel is much larger than the experimental value and the 
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bound state is observed in the ^Sq channel. The 2+1 flavor results from Ref. 
0.39 GeV give non-zero negative values for AEl in both channels on the <(3.9 fm)'^ box, 
which are consistent with our results as shown in Table [Vj Unfortunately, the extrapolation 
to the inflnite spatial volume limit introduces large errors so that Ai^oo becomes consistent 
with zero within the error bars. The most recent study ^ worked at a heavier quark mass 
of m^r = 0.81 GeV in 3-flavor QCD, and found large values for the binding energies: 25(3) (2) 
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MeV for the ^Si channel and 19(3)(1) MeV for the ^Sq channel ^. While all recent studies 
are consistent with a bound ground state for both ^Si and ^Sq channels when quark masses 
are heavy, quantitative details still need to be clarified. 

IV. CONCLUSION AND DISCUSSION 

We have calculated the binding energies for the helium nuclei, the deuteron and the 
dineutron in 2+1 flavor QCD with m.^ = 0.51 GeV and niN = 1.32 GeV. The bound states 
are distinguished from the attractive scattering states by investigating the spatial volume 
dependence of the energy shift AE^. In the infinite spatial volume limit we obtain 

43(12) (8) MeV for ^He, 
20.3(4.0) (2.0) MeV for ^He, 
11.5(1.1)(0.6) MeV for ^Si, 
7.4(1. 3)(0.6) MeV for ^Sq. 

While the binding energy for the ^He nucleus is comparable with the experimental value, 
those for the ^He nucleus and the deuteron are much larger than the experimental ones. 
Furthermore we detect the bound state in the ^Sq channel as in the previous study with 
quenched QCD, which is not observed in nature. These findings and the enhanced binding 
energies at m^r = 0.81 GeV in 3-fiavor QCD Q] tell us that a next step of primary importance 
is to reduce the up-down quark mass toward the physical values. A possible scenario in 
the two-nucleon channels is as follows. The binding energy in both channels diminishes 
monotonically as the up-down quark mass decreases. At some point of the up-down quark 
mass the binding energy in the ^Sq channel vanishes and the bound state evaporates into 
the attractive scattering state, while the binding energy in the ^Si channel remains finite 
up to the physical point. This is a dynamical question on the strong interaction, and only 
lattice QCD could answer it. 
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TABLE I: Simulation parameters for gauge configuration generation at {Kud, 
Ks)=(0-1373316, 0.1367526). The definition of parameters is same as in Ref. [l^. 





■^93 x 48 


40'^ X 48 


4r3 V 4Q 


fi43 X 64 


^ run 


1 


1 


1 


2 


{No,Ni,N2) 


(2,2,10) 


(2,2,15) 


(2,2,16) 


(2,2,18) 


Block size 


8^ X 6 


10^ X 6 


122 X g2 


8^ X 4 


iVpoly 


260 


320 


320 


340 


MD time 


4000 


2000 


2000 


(1090,810) 


Pacc(HMC) 


0.840 


0.925 


0.916 


(0.880,0.867) 


Pacc(GMP) 


0.957 


0.969 


0.963 


(0.978,0.974) 



TABLE II: Number of configurations, separation of trajectories between each measurement, bin 
size in jackknife analysis, number of measurements on each configuration, exponential smearing 
parameter set {A,B) in Eq. (jl2p . pion mass ttItt and nucleon mass mjy are summarized for each 
lattice size. 



L 


T 


# config. 


''"sep 


bin size 


^ meas. 


{A,B) 


[GeV] 


mN [GeV] 


32 


48 


200 


20 


10 


192 


(1.0,0.18) 


0.5109(16) 


1.318(4) 


40 


48 


200 


10 


20 


192 


(0.8,0.22) 


0.5095(8) 


1.314(4) 


48 


48 


200 


10 


20 


192 


(0.8,0.23) 


0.5117(9) 


1.320(3) 


64 


64 


190 


10 


19 


256 


(0.8,0.23) 


0.5119(4) 


1.318(2) 
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TABLE III: Energy shift —AEi in physical units and fit range for ^He and ^He channels on each 
spatial volume. Extrapolated results in the infinite spatial volume limit are also presented. The 
first and second errors are statistical and systematic, respectively. 





^He 




^He 




L 


-AEl [MeV] 


fit range 


-AEl [MeV] 


fit range 


32 


47(24) (5) 


10-14 


23.2(7.6)(1.4) 


10-14 


40 


30(15)(23) 


9-13 


20.2(6.9)(2.8) 


9-14 


48 


39(20)(27) 


10-14 


25.5(5.3)(1.7) 


10-14 


64 


46(11)(8) 


10-14 


19.5(3.7)(1.2) 


9-14 


oo 


43(12)(8) 




20.3(4.0)(2.0) 





TABLE IV: Same as Table Hill for ^Si and ^Sq channels. 



L 


-AEl [MeV] 


fit range 


-AEl [MeV] 


fit range 


32 


12.4(2.1)(0.5) 


9-14 


6.2(2.4)(0.5) 


10-14 


40 


12.2(1.9)(0.6) 


10-15 


8.2(4.0)(1.5) 


11-15 


48 


11.1(1.7)(0.3) 


10-14 


7.3(1. 7)(0.5) 


10-14 


64 


11.7(1.2)(0.5) 


9-14 


7.2(1.4)(0.3) 


10-14 


oo 


11.5(1.1)(0.6) 




7.4(1.3)(0.6) 
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TABLE V: Energy shift —AEl in physical units for ^Si and ^Sq channels together with the previous 
works. The values marked by * are estimated from the scattering length in Ref. Q] employing the 
leading term of finite volume formula in the 1/L expansion [2^ with the nucleon mass obtained 



from the same ensemble in Ref. 
interpolating operator. 



30| . The values for Ref. are taken from the results with the Oi 



Ref. quark action # flavor [GeV] L [fm] -AEl [MeV] 











01 




[71 Wiknn 





n 79 


9 7 

Li. I 


99 ^\ 

Zic/ . 1 U. t? j 


1 4 7(A '^'i 




n 
u 


D QQ 


9 7 


LO . 1 I U.O 1 


1 7(A 

±U. ( t '4:.0 I 







1 


9 7 

L. 1 


XO . ± I U . U 1 


1 9 9/"^ Q^l 

± Lj . Lj\0 .cf J 


iViiJvcCL [LJ VV Uii YT-oV^tdU. 1 


9-1-1 


\j .00 




— 1 fiM Q^* 


XUl -LO ) 




9-1-1 


n zLQ 




— y .0 i^u.tj j 


1 ^ 1 l^zL 9^* 




9-1-1 


n 

U. eJt? 


9 




u.ui 1 . I 


[111 Wilt^nn 


n 


u.oo 


d A 


U. t7 M 1 ) 


U.UOl ZrU 1 




n 


U. (JO 


A A 




U.UUc7l^e7'± j 




n 


u. / 


A A 


u.'4:OUi^y 1 ^ 


U.'4:UUl^OO ^ 


V \ A/ 1 1 c r\i~i c cw Tc^y 
VV iloUll-L^iU Vci 


n 


n Ro 


1 


1 n 91^9 9V1 fi\ 


U.llisi.O }\Li.Li J 




n 


u.ou 


1 

U. 1 








n 


D RO 
u.ou 


19*^ 

LL.O 


t .0{Li. 1 I I U.^^: ) 


A f<(o nVi 1 




n 


n Ro 

u.ou 




y. 1 1 1 . 1 J I u.u ) 


tJ.tJlJ-.lMX.U) 


1 J 1 jTiiilbO. VV lolUll-CiUVci 


9-1-1 


u.oy 


9 A 


1 f^{9 f^\(A '^^ 
l.yjyLi.K) J y^.O J 


Q^^ 7\(o fi'^ 

0. yi^i . / j yLi.xj j 


[3] Aniso. Wilson-clover 


2+1 


0.39 


3.0 


22.3(2.3)(5.4) 


10.4(2.6)(3.1) 




2+1 


0.39 


3.9 


14.9(2.3)(5.8) 


8.3(2.2)(3.3) 




2+1 


0.39 


00 


11(5)(12) 


7.1(5.2)(7.3) 


[4] Stout Wilson-clover 


2+1 


0.81 


3.4 


25(3)(2) 


16(3)(1) 




2+1 


0.81 


4.5 


21(3)(1) 


11(2)(1) 




2+1 


0.81 


6.7 


25(3)(2) 


19(3)(1) 


This work Wilson-clover 


2+1 


0.51 


2.9 


12.4(2.1)(0.5) 


6.2(2.4)(0.5) 




2+1 


0.51 


3.6 


12.2(1.9)(0.6) 


8.2(4.0)(1.5) 




2+1 


0.51 


4.3 


11.1(1.7)(0.3) 


7.3(1.7)(0.5) 




2+1 


0.51 


5.8 


11.7(1.2)(0.5) 


7.2(1.4)(0.3) 




2+1 16 


0.51 


00 


11.5(1. 1)(0.6) 


7.4(1.3)(0.6) 



0.61 



0.605 - 



0.6 - 



0.595 - 



m 



N 



12 



16 



20 



FIG. 1: Nucleon effective mass on (5.8 fm)^ box in lattice unites. Fit result with one standard 
deviation error band is expressed by solid lines. 




FIG. 2: Effective energy shift AEJ^ for ^He channel on (5.8 fm)^ box in lattice units. Fit result 
with one standard deviation error band is expressed by solid lines. 
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FIG. 3: Spatial volume dependence of AEl in GeV units for ^He channel. Outer bar denotes 
the combined error of statistical and systematic ones added in quadrature. Inner bar is for the 
statistical error. Extrapolated result in the infinite spatial volume limit is shown by filled square 
symbol together with the fit line (dashed). Experimental value (star) and quenched result (open 
diamond) are also presented. 




FIG. 4: Same as Fig. [2] for ^He channel. 



18 



-0.02 



-0.04 - 



-0.06 - 



-0.08 



AE^( He) [GeV] 



■Jl^ experiment 
9 Nj=2-hl m^=0.51GeV 
N =0 m =0.80GeV 

>^ f n 



le-05 



2e-05 



3e-05 



4e-05 



FIG. 5: Same as Fig. [3] for ^He channel. 
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FIG. 6: Same as Fig. [2] for '^Si channel. 
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FIG. 7: Same as Fig. [2] for ^Sq channel. 
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FIG. 8: Same as Fig. [3] for '^Si channel. 
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FIG. 9: Same as Fig. [3] for ^Sq channel. 
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FIG. 10: dependence of AEo^ for ^Si channe 

2+1/3 flavor (quenched) result. The results of Refs. 



.. Closed(open and cross) symbol denote the 
2, l3| and this work are extrapolated values in 



the infinite volume limit. Experimental result (star) is also presented for comparison. 
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FIG. 11: Same as Fig. [TOlfor ^Sq channel. 
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